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We describe parton density function sensitivity studies, using muon final states produced through 
the Drell-Yan process via W, Z or y* down to a Q 2 of 10 GeV 2 . This makes use of LHCb's unique 
Q_i ability to trigger on low transverse momentum objects. Due to the forward acceptance of LHCb, 

x values down to 2xl0~ 6 can be probed, where with just 100 pb^ 1 of data, the gluon PDF can 
improved up to 70%. 
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1. Introduction 

LHCb [1] is a detector designed to study b-physics so its pseudo-rapidity coverage extends in 
the forward region where B mesons are dominantly produced (1.9 < f] < 4.9). Part of that range 
is in common with ATLAS and CMS which have a coverage of |t]| < 2.5; the rest is unique to 
LHCb. During the first year of data taking, LHC will run at the center of mass energy of 7 TeV 
and an integrated luminosity of 1 fb~ l will be delivered to LHCb in the first year of data taking. 
A unique feature of LHCb is its ability to trigger and reconstruct very low transverse momentum 
muons (pj > IGeV). Besides b-physics goals LHCb is able to perform precise measurement of 
electroweak channels with muon final states at high rapidities. We will show that these measure- 
ments can probe parton distribution functions for different values of Bjorken-x. 



2. LHCb parton kinematics 



The current knowledge of the PDFs arises from fit- 
ting procedures to HERA , fixed target experiments and 
Tevatron data (or a subset of those). The foreseen elec- 
troweak measurements at LHCb include the analysis of 
the Z and W ± production and the Drell-Yan pair produc- 
tion of muons in four different mass bins (2.5 — 5 GeV , 
5-10 GeV, 10 - 20 GeV and 20 - 40 GeV )[2]. Pro- 
duction at high rapidity involves the interaction of one 
valence quark carrying a large fraction of proton mo- 
mentum (high-x) and a sea quark with small fraction of 
momentum (very small x); according to the relation that 
link x and Q 2 
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Figure 1: Percentage uncertainty on pro- 
duction cross-section of electroweak bo- 
son as function of rapidity 



where -Js is the energy in the center of mass of the colli- 
sions. For the Drell-Yan production with Q in the region 
2.5 — 5 GeV the lowest value of x reachable is ~ 2xl0~ 6 . 
Hence due to its forward rapidity coverage and the abil- 
ity to reconstruct low invariant mass muon pair systems, LHCb will cover an un-explored region of 
kinematical space where Parton Distribution Functions are known with less precision. This aspect 
can be also deduced looking at the theoretical uncertainties on the rapidity distribution for W , and 
Z/y*; The largest contribution to these uncertainties arises from PDFs uncertainties. In figure 1 
the percentage uncertainties on the production cross-sections for W , Z and Drell-Yan process at 
three different mass are reported: the uncertainties due to PDF become large at high rapidity for 
all the channels studied (corresponding to rapidity LHCb acceptance) and for the Drell-Yan pair 
production, the smaller the mass, the larger the uncertainty. 



3. The fitting algorithm 

PDFs are obtained via global fits to experimental data by different working groups, and their 
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results have some model-dependence that may differ between each other. Theoretical predictions 
for the PDF sets MSTW08[3], CTEQ66[4], Alekhin[6] and NNPDF2.0[5] are used in this analysis. 
The first three groups perform a global fit based on a Hessian method while the fourth uses a neural 
network approach based on Montecarlo replicas. We now briefly describe our method to constrain 
the PDFs where the PDFs set have been obtained using these two approaches. 

3.1 PDFs built using a Hessian method 

The Hessian method is used by the groups MSTW, CTEQ and Alekhin. A polynomial parametriza- 
tion of the PDFs such as xf(x, <2o) = (1 — -^(l + £* ' 5 + yx)x s is chosen. The global minimum 
of their fits to the global data produces a central value with the uncertainties given by a N x N 
covariance matrix. Diagonalizing the matrix gives a set of N independent variables. In this way a 
representation of the PDF in a space where uncertainty can be treated as independent with a set of 
N eigenvectors is provided. 

A given physical observable, for example the rapidity distribution, f(y), can be described as 
a linear combination of the observable obtained using the PDF best-fit, fo(y), and the observable 
fi(y) obtained by moving the i ,h eigenvector by 1-sigma. 

f(y) = My) + £ [My) - fi(y)} A,- (3. l) 

i=i 

The one sigma bound is mapped out by sampling A,- from a unit Gaussian with mean zero; this 
function also provides a form with which the data can be fitted. Given binned data with Nj events 
in the / bin with error Sj, the normalization and the shape of the distribution can be fitted by 
minimizing 



#bins 
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where Ao is the overall normalization or luminosity and A,- are the fitted eigenvalues, constrained 
by the second term to follow a Gaussian probability. Before the fit the eigenvalues by definition of 
the model have an uncertainty equal to 1. After our fit, values smaller than 1 imply that we have 
a better constraint on the PDF. The hessian method provides a diagonalized covariance matrix of 
the parameters. This fitting operation provides a new covariance matrix, V, in which values on the 
main diagonal are smaller than 1 and terms off diagonal can be different from zero. This matrix 
can be used to calculate uncertainties on PDFs using the usual equation for correlated variables: 



/_ dF dF 

where A, are the fitted parameters and F is a PDF (e.g.e uy, dy , sea , gluon etc.). 

3.2 PDFs built using a Neural network method 

NNPDF uses a neural network approach. PDFs are parametrized with a redundant flexible 
neural-net (i.e. a convenient functional form with ~ 200 parameters) This technique produces a set 
of equally likely Montecarlo replicas of the PDFs. The best estimation of any observable of the 
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PDFs is obtained by generating this observable using each of the replicas and averaging: 
(f[l]) = j/Hk=if [l 1 where N is the number of replicas, / is a generic function of the PDF 
set and qW are PDF replicas. The uncertainty can be expressed through the standard deviation. 
°/ = (w=T £/t=i (/ W®\ ~ (f M)) 2 ) ■ ^ more than one observable is generated, it is even possible 
to compute correlations between quantities. 

To evaluate the PDF constraint with this PDFset we generated pseudo-experiment as average over 
the replicas then we generated rapidity distributions with each replica. We compared the prediction 
of a single replica with the pseudo-experiment minimizing the x 2 

2 _y f Nbin-hofi(y) 

%i A,. 

bin V ^ bm 

where i is the index of the replica, Nbm is the pseudo-experiments result in that rapidity-bin and fj(y) 
is the value of the rapidity distribution generated with the i th replica. This algorithm associates to 
each replica of the PDFs a % 2 and so the degree of compatibility of each replica can be evaluated as 
the chi-square probability. Now it is possible to constrain the PDFs considering the ^-probability 
as the conditional probability for measuring the set of NNPDF replicas considering the new LHCb 
data [7]. We don't consider all the PDF replicas as equally likely anymore but we use the % - 
probability as a weight in the central value and uncertainties calculations. 



4. Results and comments 

We present the expected improvement on u-valence, d-valence, sea , s + and gluon distributions 
when introducing 1 fb~ l of LHCb electroweak data. Here we present results obtained using the 
PDF set MSTW08. Results for the other PDF sets can be found in [8]. In figure 2 the ratio 
between the PDF uncertainty after and before the fitting procedure is plotted. As expected the 
largest improvement is given by the lowest mass bin of the Drell-Yan process. With 1 fb~ l of 
LHCb data the sea and the gluon distributions at x ~ 10~ 5 and at high-x can be improved up to 
70% while in mid-x regions the improvement is about 50%. Other PDFs can be improved on the 
full x range of about 60%. For the W ± and Z production the best result is obtained when performing 
a combined fit of the three rapidity distributions. About 30% improvement can be expected over 
the whole x range for each of the PDFs. 
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(a) Improvement when fitting W + , W Z rapidity distributions and a combined fit 
of the three 
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(b) Improvement when fitting Drell-Yan rapidity distributions in 4 different mass 
bins 

Figure 2: Ratio between PDF uncertainties before and after adding 1 fb~ 1 of LHCb data for different 
distributions as a function of x 
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